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Abstract. In order to analyze spectrum from the interstellar medium (ISM), spectrum of 
the molecule of interest is recorded in a laboratory, and accurate rotational and centrifugal 
distortion constants are derived. By using these constants, one can calculate accurate 
partition function. However, in the same paper, where these constants are derived, the 
partition function is calculated by using a semi-empirical expression. We have looked into 
the details of this semi-empirical expression and compared the values, obtained from it, 
with the accurate ones. As an example, we have considered the case of Methanimine 
(GH 2 NH) which is detected in a number of cosmic objects. It is found that for the kinetic 
temperature T > 120 K, the semi-empirical expression gives large value as compared to 
the accurate one. The deviation becomes about 25% larger than the accurate one at the 
kinetic temperature of 400 K. 
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1 Introduction 

Partition function derived from the spectroscopic parameters connects to the thermody¬ 
namic quantities. From time to time, in various contexts in Astronomy & Astrophysics, 
computation of partition function has been a point of extensive discussion. For example, 
the calculation of partition function of Hg" has been in hot discussion in the context of 
solar continuum observations. Several scientists calculated the partition function of Hg". 
The partition function of Hg" as a function kinetic temperature, calculated by Chandra 
et al. (1991), was declared by Sidhu et al. (1992) to be the best one. Later on, Neale & 
Tennyson (1995) extended the calculation to high temperatures. The rotational partition 
function Q is expressed as (Chandra et ai, 1991; Sidhu et al., 1992; Neale &: Tennyson, 
1995) 

Q = Y.{2Ji + 1) eM-hcEi/kT) (1) 
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Here, T is the kinetic temperature and the nuclear degeneracy is not considered. The h, 
k, c are, respectively, the Planck constant, Boltzmann constant and speed of light. For 
a rotational level i of asymmetric top molecule, the parameters Jj, Ei are the rotational 
quantum number and energy, respectively. The (2Jj + l) is the rotational statistical weight 
of level i. The values of Ji and Ei for a level can be calculated by using the rotational and 
centrifugal distortion constants for the molecule. 

However, in many papers where rotational and centrifugal distortion constants are 
derived with big efforts, the partition function is calculated by using the semi-empirical 
expression 



( 2 ) 


where A, B and C are rotational constants. Details of this expression (2) are discussed 
in the next section. In the present paper, we have discussed about the deviation obtained 
when one uses equation (2) in place of (1). 


2 About equation (2) 


Kassel (1933) used partition function for calculation of thermodynamic quantities. Asymp¬ 
totic expansion of the expression for the rotational partition function of a rigid symmetric 
top molecule for high temperature T was given by Viney (1933) as the following. 
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where A and B are rotational constants for the symmetric top molecule. 

Since there is no explicit formula for energy of rotational level in case of an asymmetric 
top molecule, it is impossible to derive a rigorous asymptotic expansion for the partition 
function in this case. However, Gordon (1934) suggested that for a rigid asymmetric top 
molecule, in equation (3), B should be replaced by VBC to get the expression as 


Q = 
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Here, A, B and C are rotational constants for the asymmetric top molecule. This equation 
(4) may be considered as a semi-empirical expression for a rigid asymmetric top molecule. 
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For sufficiently high temperature (or small rotational constants), equation (4) reduces 
to (2), as the argument of the exponential would be very small, and the contributions 
of second, third, etc. terms would be neglected. This expression (2) is thus a semi- 
empirical relation derived for a rigid asymmetric top molecule and is valid at sufficiently 
high temperature. Moreover, the expression (4) is better than (2). But, this expression 
(2) is being used regularly by the scientists. Probably they are not aware of the facts 
about the expression (2), and as a tradition, this expression is being used. 

3 Results and discussion 

Now, we have three equations (1), (2) and (4), for the partition function. Let us denote 
their values as Qi, Q 2 , Q 3 , respectively, for convenience. There are many papers where 
accurate rotational and centrifugal distortion constants of asymmetric top molecule are 
derived, and in the same paper the partition function is calculated with the help of equation 
(2). Some examples of such papers are: Widicus et al. (2003) (for NH 2 CH 2 CH 2 OH), 
Widicus Weaver et al. (2005) (for CHOCH 2 OH), Motoki et al. (2014) (for CH 2 NH). We 
have not intention about them. We have considered the case of Methanimine (CH 2 NH) 
and have used the data of Motoki et al. (2014), given in Table 1. 

Using the values given in Table 1, we have calculated accurate energies of rotational 
levels with the help of the computer code ASROT (Kisiel, 2001). For calculation of 
partition function, we have considered the levels having energy up to 1475 cm“^, which 
is sufficiently large for the kinetic temperature up to 400 K, considered here. Using these 
energies of levels, the partition function is calculated with the help of equation (1). At 
400 K, the contribution of the highest level considered is less than 0.01%. It shows that 
the considered set of energy levels is sufficient, and further high energy levels would not 
contribute to the partition function significantly, up to a temperature of 400 K. 

The results obtained Qi and Q 2 are presented in Figure 1 as a function of kinetic 
temperature T in the gas. The values of partition function reported by Motoki et al. (2014) 
were obtained after multiplying by a factor of 3 for the nuclear degeneracy. Therefore, 
they differ by a factor of 3 as compared to those obtained with the help of equation (2). 

As obvious, the Q 3 is always larger than Q 2 . For the temperature T < 120 K, the 
deviation of Q 2 with respect to Qi is less than 1%, and Q 2 is always smaller than Qi. 
For T > 120 K, the deviation of Q 2 with respect to Qi increases with temperature. The 
deviation becomes 25% at T = 400 K. It is quite large, and the values of thermodynamic 
quantities would change significantly. 
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The contribution of centrifugal distortion constants to partition function increases 
when we go to the levels with large value of J. In some molecules, the contribution even 
may change the sequence of levels (Sharma et al., 2014). It is interesting to note that 
equation (2) does not give accurate results even at high temperatures. 

4 Conclusion 

It may finally be concluded that when the rotational and centrifugal distortion constants 
are available, the partition function should be calculated with the help of equation (1). It 
may require some efforts, but the reliable results would be available in literature. 
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Table 1. Rotational and centrifugal distortion constants in MHz. 


Parameter 

Value 

Parameter 

Value 

A 

196210.87629 

Hjk 

2.6068 X 10-6 

B 

34641.699423 

Hkj 

6.125 X 10-6 

C 

29351.496787 

Hk 

877.75 X 10-6 

Dj 

56.00411 X 10"3 

hi 

33.067 X 10-9 

Djk 

598.6775 x lO'^ 

h2 

32.462 X 10-9 

Dk 

6.383483 

h3 

9.091 X 10-9 

di 

-9.973179 X 10-3 

Lr 

-163.0 X 10-9 

d2 

-2.007905 X 10-3 

h 

-0.541 X 10-32 

Hj 

0.01736 X 10-6 

h 

-0.472 X 10-32 


Figure 1: Variation of partition function Q versus kinetic temperature T. Solid line is for Qi, 
obtained from equation (1) and the dotted line for Q 2 , obtained from equation (2). 
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